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Domain Wall Fermions and MC Simulations of Vector Theories. 

p. Vranas^ 

''Dept. of Physics, Columbia University, New York, NY 10027 

It is known that domain wall fermions may be used in MC simulations of vector theories. The practicality 
and usefulness of such an implementation is investigated in the context of the vector Schwinger model, on a 2+1 
dimensional lattice. Preliminary results of a Hybrid Monte Carlo simulation are presented. 
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1. Introduction 

Domain Wall Fermions (DWF) for lattice 
gauge theories were developed in |^ . Since then a 
wealth of activity has followed with the main fo- 
cus being the application of DWF to chiral lattice 
gauge theories, but also to vector lattice gauge 
theories (see Q and references therein). As a 
result of these works little doubt remains that 
DWF can be used to regularize vector gauge the- 
ories on the lattice. Because DWF have many 
attractive features it is surprising that little work 
has been done in the past few years investigat- 
ing their practicality and usefulness in numerical 
simulations of lattice vector gauge theories. A 
preliminary investigation of this, in the context 
of the vector Schwinger model, is presented here. 

2. The model 

The massless case of the one flavor model is 
described by the lattice Euclidean action: 
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In the above U is the U(l) gauge field. Up the 
standard plaquette, (3 = g^^ with go the gauge 



coupling, and n is a collective space time coor- 
dinate. The photon mass is = ffo/v^ S'Ud. 
the space time volume \s V = . The key in- 
gredient of DWF is the introduction of the extra 
direction s, with periodic boundary conditions, 
size Lg, and a mass defect at s = and s = Ls/2 
{m{s) = -fmO for < s < Ls/2 and m(s) = — mO 
for Lg/2 < s < Ls,0 < mo < 1). In the free case 
a Dirac fermion appears with its right component 
exponentially peaked at s = and its left one at 
s = Ls/2. The overlap of the left and right com- 
ponents is small and decreases exponentially with 
increasing Ls- In the limit — > cx3 one obtains 
a single massless Dirac fermion. The introduc- 
tion of the extra direction also introduces heavy 
modes. Their contribution needs to be subtracted 
1^. This is done by dividing the fcrmionic deter- 
minant, det[A^] by ^det[Al+] dct[A^_], where 
M + , Ai- are the same as Ai but with m{s) — 
+mo, m(s) = —mo respectively. These determi- 
nants can be produced by introducing appropri- 
ate auxiliary bosonic fields. 

There are two implementations of DWF that 
can be used in numerical simulations: 
I The overlap formalism H). A transfer matrix 
is constructed along the extra direction and from 
it the corresponding Hamiltonian Ti. is extracted. 
This formulation allows for the strict Ls oo 
limit to be taken. Observablcs can be calculated 
by obtaining all the eigenvalues and eigenvectors 
of (Ti is a matrix of size ^ V xV). The lattice 
vector Schwinger model was simulated success- 
fully using this method Q. However, the method 
requires large amounts of computer time and high 
statistics simulations in four dimensions and large 
volumes may not be possible at present. 
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II A direct simulation of (|l|) with Lg finite. The 
chiral limit is obtained by extrapolating to the 
Lg ^ oo limit. This method is very attractive for 
several reasons: 

a) Unlike Wilson fermions there is no fine tuning 
involved; the chiral limit is the Lg ^ oo limit. 
At the Ls — > oo limit the quark mass is multi- 
plicatively renormalized |^ (for a recent investi- 
gation of this using Hybrid Monte Carlo (HMC) 
see Q). Although for finite Ls there will always 
be an additive component its size is expected to 
be decreasing exponentially with Lg. 

b) Unlike staggered fermions there is no breaking 
of flavor symmetry on the lattice. 

c) Because the overlap of the left and right com- 
ponents is exponentially small in the free case, 
one would expect that very small explicit break- 
ing of chiral symmetry can be obtained for a rel- 
atively small Lg- If this is the case, then the 
method will be very valuable in studies involving 
spontaneous chiral symmetry breaking. 

d) Anomalous symmetry breaking has a natural 
interpretation jl] and one may expect that this 
will facilitate related studies. 

Of course, there is a price to be paid for these 
nice features. The theory has one more dimension 
and therefore demands more computer resources. 

Given the above considerations several ques- 
tions should be answered before DWF can be 
used to simulate lattice vector gauge theories and 
in particular lattice QCD. Are there any hidden 
difficulties in an HMC simulation of the model? 
Does the method work and give the correct an- 
swers? How does the computational difficulty de- 
pend on Lg and how large should Lg be? How 
does the computational difficulty of DWF com- 
pare with Wilson or Staggered fermions? How 
good are DWF in addressing questions related 
to anomalous breaking of axial symmetries and 
spontaneous breaking of chiral symmetries? Pre- 
liminary answers to a few of these questions are 
presented below. 

3. The size of the extra dimension 

The Lg dependence can be investigated by di- 
rect comparison with the Lg oo result obtained 
using the overlap formalism. In the massless case 



the fermionic effective action Sgf / is finite in the 
zero topological sector and Seff — *■ — oo other- 
wise. On an 8 X 8 lattice, Seff is calculated for 
several Lg for a background gauge field with topo- 
logical charge in figure 1, and with charge 1 in 
figure 2. The very different Lg dependence in the 
two topological sectors is evident. This indicates 
that DWF correctly reproduce topological effects 
already at Lg ~ 10—12. Furthermore, the overlap 
result in the zero topological sector (dashed line 
in figure 1) is also reached at Lg = 10 — 12. Sim- 
ilar results for the size of Lg have been obtained 
from studies of the pion mass in the sector B. 




Figure 1. Seff in topological sector 
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Figure 2. Seff in topological sector 1 



4. Hybrid Monte Carlo and anomalous 
symmetry breaking 

A variant of (|l]) developed in |^ is most suited 
for use with the HMC algorithm. The basic 
difference is that m(s) — niQ but free bound- 
ary conditions are used along the extra direction. 
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The two chiral modes are now bound to the two 
free boundaries. A mass term connects the free 
boundaries and introduces a bare quark mass nif. 
The advantage of this is that it only needs half 
the size of Lg and also only a single bosonic field 
to perform the subtraction mentioned earlier. 

The massive Nj = 2 flavor vector Schwinger 
model is simulated using this formulation and 
the HMC algorithm. The expectation value of 

the operator W = Uf'Ji *k*L + Il^I^i *L*fl is 
calculated. If the volume is kept fixed and mf 
is made very small the effect of the zero modes 
coming from different topological sectors becomes 
important 0. In particular, as can be seen from 
the overlap implementation Q [Q , in the massless 
and Ls oo limit the operator W receives con- 
tributions only from sectors ±1, while the fermion 
Boltzman weight (fermion determinant) is not 
zero only in sector 0. As is turned on (and/or 
Ls is decreased from infinity) the fermionic deter- 
minant becomes non zero in sectors other than 
and the operator W receives contributions from 
sectors other than ±1. Therefore for small m/ 
the HMC algorithm will mostly sample the sec- 
tor where the observable W receives small con- 
tributions. The algorithm will infrequently visit 
the sectors ±1 but when it does the observable 
W will receive large contributions to make up for 
the small sampling rate. As a result, when my is 
decreased a larger number of HMC iterations will 
be needed to sample the ±1 sectors correctly. 
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HMC sweep 
Figure 3. Time history oi W a^t mf — 0.01 



A 
V 



0.05 



0.04 



0.03 



o = 6 
□ = 8 
X overlap 



0.0 



0.2 



0.4 



Figure 4. (W) vs. m/ 
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5. Preliminary results 

Preliminary results are presented for mg = 0.9, 
L = 6, m^L — 3.0. The HMC sweeps vary from 
6,000 for m/ = 0.01 to 2,000 for m/ = 0.5. The 
trajectory length is 1 and the step size is 0.02 for 
mf — 0.01 and 0.04 for the rest. The average 
conjugate gradient iterations vary from 56 to 86. 
The time history oi W for mf — 0.01 is shown 
in fig. 3. The large "spikes" are related to con- 
figurations with charge ±1 (measured with the 
geometric method). (W) vs. mf is in fig. 4 and 
agrees roughly with the m/ = overlap result [|j . 
For comparison with the exact answer see ^ . 
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